%This function calculates the junction area as a function of the width of
%the tape, w, and the junction angle, phi
function [Areaslist]=areacalculation(w,testangles)
Areaslist=zeros(1,length(testangles));
for i=1:length(testangles)
    testangle=testangles(i);
    Area=@(phi) (1./2).*(w.^2).*((secd(phi./2)).^2).*(cotd(phi./2));
    Areaslist(i)=Area(testangle);
end